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摘 ES 针对 全 球 定位 系统 (Global Positioning System, GPS) 星 载 原子 钟 在 钟 差 预 报时 


与 不 同 模 型 的 适应 度 不 同 的 问题 ， 采 用 二 次 多 项 式 CQuadratic Polynomial, QP) 模型 、 灰 色 

(Grey Model, GM(1,1)) 模型 和 灰色 + 自 回归 ( GM(1,1)+ Autoregressive, AR) 模型 对 不 同类 型 
原子 钟 的 钟 差 进行 预报 ,着 重 分 析 不 同类 型 原子 钟 的 预报 精度 、 不 同 长 度 钟 差 序 列 建 模 预报 
效果 以 及 钟 差 序 列 波动 对 预报 结果 的 影响 。 实 验 结果 表明 : 钟 差 预报 精度 与 建 模 序列 长 度 有 
一 定 关 系 ， 二 次 多 项 式 模型 受 影响 最 大 ， 灰色 + 自 回 归 模 型 受 影响 最 小 ， 不同 卫星 原子 钟 在 
不 同 预报 模型 下 最 佳 建 模 序 列 长 度 不 同 , 鲍 钟 受 建 模 序列 长 度 的 影响 小 于 钨 钟 ; 二 次 多 项 式 
模型 对 钨 钟 预报 效果 较 差 ， 对 鲍 钟 预报 效果 可 与 灰色 模型 和 灰色 + 自 回 归 模 型 相当 ; 钟 差 序 
列 波动 时 ， 建 模 预报 精度 降低 ， 不同 模型 的 预报 结果 受 钟 差 波 动 幅度 大 小 的 影响 不 同 。 
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卫星 钟 差 (Satellite clock bias, SCB) 是 全 球 卫 星 导 航 系 统 (Global Navigation Satellite 
System, GNSS) 实现 高 精度 定位 、 导 航 和 授时 的 重要 参数 之 一 ， 但 国际 GNSS 服务 
(International GNSS Service, IGS) 发 布 的 精密 钟 差 产品 会 有 一 定 程度 的 延迟 上 3， 为 实现 用 
户 实时 、 高 精度 的 使 用 要 求 ， 对 卫星 钟 差 的 预报 进行 研究 就 成 为 当前 的 重要 任务 之 一 。 

用 于 卫星 钟 差 预 报 的 模型 主要 有 二 次 多 项 式 模型 中、 灰色 模型 个、 自 回 归 滑 动 平均 
(Autoregressive Moving Average, ARMA) AIO, Kalman 滤波 中 和 人 工 神经 网 络 模 型 9 
等 ,不 同 模型 在 预报 时 均 有 各 自 的 优势 和 局 限 性 外 。 由 于 全 球 定位 系统 星 载 原子 钟 的 类 型 不 
完全 相同 ， 其 特性 及 稳定 性 也 有 所 不 同 ， 故 其 适用 的 预报 模型 也 不 一 定 相 同 中 。 不 少 学 者 
对 钟 差 进行 预报 时 , 主要 是 验证 模型 的 预报 精度 , 并 没有 对 不 同类 型 原子 钟 更 适合 哪 种 模型 
进行 分 析 。 另 外 ， 在 研究 预报 模型 的 过 程 中 ， 多 是 选取 固定 长 度 的 原始 钟 差 序列 进行 预报 ， 
鲜 有 考虑 采用 不 同 长 度 的 原始 钟 差 序 列 建 模 对 预报 精度 的 影响 , 且 多 数 实验 均 是 选择 变化 趋 
势 较 为 稳定 的 卫星 钟 差 数据 , 也 没有 对 钟 差 具 有 一 定 波动 情况 下 不 同 模 型 的 预报 精度 进行 分 
析 。 


本 文选 取 灰 色 模 型 、 灰 色 + 自 回归 模型 和 二 次 多 项 式 模型 对 全 球 定位 系统 卫星 钟 差 的 短 
期 预报 进行 研究 ， 首 先 ， 分 析 采 取 不 同 预报 模型 时 ， 不 同类 型 原子 钟 的 预报 精度 ， 其 次 ， 对 
不 同类 型 原子 钟 预报 时 , 选取 不 同 长 度 的 钟 差 序列 进行 建 模 ， 研 究 不 同 卫星 在 不 同 模型 下 的 
最 佳 建 模 长 度 ， 最 后 ,选取 钟 差 序列 波动 的 时 间 段 进行 预报 , 分 析 钟 差 序 列 波动 大 小 对 不 同 
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模型 预报 精度 的 影响 。 
1 预报 模型 
1.1 二 次 多 项 式 模型 
若 原始 钟 差 序 列 为 X = (人 x 0x) ， 则 二 次 多 项 式 模型 可 表示 为 : 
X =a +a (t; —t)-*aj,(t, 1) +E (D 


其 中 ， 二 1,2,...,n; to 为 参考 时 刻 ; ;为 历 元 时 刻 ; ao、al 和 a 分 别 为 时刻 的 钟 差 、 钟 速 和 

钟 漂 ; ”为 误差 。 当 原始 钟 差 序列 长 度 n 大 于 3 时 ， 即 可 根据 最 小 二 乘 原理 解 算出 参数 。 

1.2 灰色 模型 
若 原 始 钟 差 序列 和 = 人 x … 二) 为 非 负 序列 , 为 增强 序列 的 规律 性 , 对 X 进行 一 次 


k 
累加 可 生成 序列 X= xe? s x0), JOB: x? -Y x 0-123..0. 
i-l 
XpxO 建立 微分 方程 GR 〈3))， 即 可 根据 最 小 二 乘 原理 即 可 求 得 参数 。 
dx? w 
tax 
dt 
其 中 ，a 为 模型 的 发 展 系数 ; b 为 灰色 作用 量 。 
依据 (2) 式 的 响应 函数 ， 可 得 X0 的 预报 模型 为 : 
b 
(1) k 


ds mx Ps 十 一 (3) 
a a 


其 中 ,大 表示 参与 建 模 的 原始 数据 个 数 ,对 预报 的 xO (Kc) 序列 依次 执行 累 减 即 可 得 到 原始 序 


列 的 预测 结果 。 
1.3 自 回 归 模 型 
若 有 平稳 时 间 序 列 X = 人 x, … 过)， 则 自 回 归 模 型 可 表示 为 品 1 


=b (2) 


E= X,-aX,,*0oX, tg X, ,+ (4) 
O HB. ee HERRA p 为 模型 阶 数 ;， e, 为 白 品 声 。 


自 回归 模型 阶 数 通常 由 AIC (Akaike information criterion, AIC? 和 BIC (Bayesian 
information criterion, BIC) 确定 ， 之 后 即 可 根据 最 小 二 乘法 解 算出 模型 系数 。 本 文采 用 BIC 
确定 模型 最 佳 阶 数 ，BIC 定义 如 下 1; 

BIC(p) = pln(n) - 21n(L) (5) 
RP, p 为 模型 阶 数 ，n 为 样本 量 ; 工 为 似 然 函 数 。 


2 实验 分 析 


为 分 析 不 同 因 素 对 钟 差 预报 精度 的 影响 , 本 文选 取 IGS 发 布 的 2016 年 10 月 01 日 至 10 
H 31 日 共 31 天 的 精密 钟 差 产品 作为 实验 数据 ， 采 样 间隔 为 Smin。 考 虑 到 全 球 定位 系统 星 
载 原子 钟 有 不 同类 型 5 主要 选取 具有 代表 性 的 PRN02 AR RB), PRN06 (IIF RB), PRN08 
(IF CS)、 和 PRN29 (CIR-M RB) 进行 实验 ， 采 用 二 次 多 项 式 模型 、 灰 色 模 型 和 灰色 + 自 
回归 模型 3 种 预报 算法 分 别 进行 预报 以 研究 原子 钟 类 型 、 原始 钟 差 序列 长 度 及 其 波动 情况 对 
预报 精度 的 影响 ， 采 用 最 大 值 (MAX)、 最 小 值 (MIN)、 平 均值 (MEAN) 和 均 方 根 误差 


Hi 


(RMS) 对 预报 结果 的 精度 进行 统计 。 
2.1 不 同类 型 原子 钟 预报 精度 分 析 
全 球 定位 系统 星 载 原子 钟 种 类 不 同 , 为 分 析 不 同类 型 原子 钟 的 预报 精度 , 分 别 以 前 一 天 
的 钟 差 (24h，288 个 历 元 〉 作为 原始 钟 差 序列 建 模 ， 预 报 接 下 来 时 长 为 6bh，12h 和 24h 的 
钟 差 , 如 采用 10 月 01 日 24h 的 钟 差 数 据 建 模 对 10 月 02 日 前 6h，12h 和 24h 的 钟 差 进 行 预 
jJ, 重复 预报 30 次 ， 并 对 30 次 预报 结果 进行 统计 分 析 。 由 于 搭载 狗 钟 的 卫星 较 少 ， 除 上 述 
4 颗 卫 星 外 ， 增 加 PRN24 OIF CS) 一 同 进 行 实验 。 由 于 篇 幅 问 题 ， 文 中 只 给 出 了 5 颗 卫 星 
24h 的 钟 差 预报 结果 ， 见 图 1 至 图 5， 表 1 给 出 了 5 颗 卫 星 不 同情 况 下 的 预报 精度 统计 
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图 3 PRN08 24h 预报 误差 图 4 PRN24 24h 预报 误差 


Fig.3 24h Prediction Errors of PRNO8 Fig.4 24h Prediction Errors of PRN24 
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图 5 PRN29 24h 预报 误差 
Fig.5 24h Prediction Errors of PRN29 
R1 3 种 模型 预报 不 同时 长 精度 统计 (单位: ns) 
Table 1 Accuracy Statistics of Three Models with Different Prediction time length (unit:ns) 
预报 模 6h 12h 24h 


型 MAX MIN MEAN RMS MAX MIN MEAN RMS MAX MIN MEAN 


PRNO2 QP 1.960 1.576 1.738 1.740 1.960 1.576 1.760 1.761 2.079 1.576 1.821 


RMS 


1.824 
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GM(1,1) 3.208 2.932 3.070 3.071 3.208 2.570 2.931 2.936 3.208 1.742 2.563 2.599 


GM(1,1) 
1.690 1.398 1.539 1.540 1.690 1.104 1.423 1.432 1.690 0.334 1.093 1.161 
+AR 
QP -1.440 -1.630 -1.531 1.532 -1.440 -1.704 -1.566 1.567 -1.440 -1.813 -1.640 1.643 
GM(1,1) -3.063 -3.211 -3.129 3.130 -3.063 -3.434 -3.221 3.223 -3.063 -3.946 -3.448 3.457 
PRN06 
GM(1,1) 
-1.501 -1.642 -1.563 1.563 -1.501 -1.891 -1.663 1.667 -1.501 -2.432 -1.905 1.924 
+AR 
QP 0.789 0.202 0.490 0.512 1.584 0.202 0.628 0.689 4.454 0.202 1.843 2.280 
GM(1,1) 0.553 -0.381 0.013 0.232 0.553 -1.447 -0.378 0.659 0.553 -1.680 -0.574 0.774 
PRNO8 
GM(1,1) 
0.669 -0.316 0.109 0.272 0.669 -1.424 -0.313 0.651 0.669 -1.679 -0.539 0.768 
+AR 
QP 0.214 -3.216 -1.736 2.013 0.214 -5.600 -3.137 3.547 0.214 -10.484 -5.400 6.099 
GM(1,1) 1.199 -1.003 -0.168 0.664 1.199 -1.645 -0.728 1.034 1.199 -1.645 -0.741 0.926 
PRN24 
GM(1,1) 
0.352 -1.732 -0.927 1.106 0.352 -2.415 -1.489 1.655 0.352 -2.443 -1.591 1.689 
+AR 
QP 0.676 -0.460 -0.099 0.348 0.676 -0.624 -0.257 0.391 0.676 -1.927 -0.839 1.069 
GM(1,1) 1.164 0.041 0.404 0.520 1.164 -0.082 0.266 0.387 1.164 -1.263 -0.262 0.659 
PRN29 
GM(1,1) 
0.488 -0.407 -0.051 . 0.258 0.488 -0.407 -0.060 0.203 0.488 -1.303 -0.458 0.653 
+AR 
通过 分 析 可 知 : 


a) 从 3 种 算法 预报 误差 的 波动 情况 看 ， 锦 钟 (PRN02，PRN06，PRN29) 误差 波动 
幅度 均 在 3ns 以 内 ， 部 分 情况 的 波动 误差 在 lns，3 种 算法 的 预报 误差 波动 幅度 类 似 ， 主 要 
于 锦 钟 较为 稳定 ; ftt (PRN08，PRN24) 的 预报 误差 波动 幅度 受 算法 影响 较 大 ， 二 次 多 
项 式 模型 的 预报 误差 波动 明显 较 大 ， 可 达 12ns， 灰 色 模 型 和 灰色 + 自 回 归 模 型 模型 的 预报 误 
差 波 动 幅 度 可 控制 在 3ns， 另 外 对 3 种 算法 预报 结果 的 平均 值 和 均 方 根 误差 对 比 ， 也 说 明 灰 
色 模 型 和 灰色 + 自 回 归 模 型 对 于 铭 钟 的 预报 相对 二 次 多 项 式 模型 更 稳定 。 

(2) 钨 钟 虽然 稳定 性 低 于 锦 钟 ， 但 从 预报 结果 的 平均 值 和 均 方 根 误 差 可 知 ， 钨 钟 采用 
灰色 模型 和 灰色 + 自 回 归 模 型 的 预报 精度 可 与 鲍 钟 相当 甚至 部 分 情况 下 由 于 锦 钟 ， 这 也 从 侧 
面 说 明了 灰色 模型 的 抗 干扰 能 力 较 强 , 而 二 次 多 项 式 模 型 预报 饮 钟 的 误差 较 大 , 说 明 二 次 多 
项 式 模型 更 适合 对 变化 趋势 较为 稳定 的 钟 差 进行 预报 。 

G) 对 比 不 同 卫星 预报 结果 ， 总 体 来 说 ， 当 原始 序列 长 度 为 24h 时 ， 采 用 灰色 + 自 回 
归 模型 的 预报 效果 最 好 , 这 是 因为 组 合 模型 从 一 定 程度 上 继承 了 单一 灰色 模型 的 优点 且 降 低 
了 灰色 模型 的 预报 风险 ， 而 AR 模型 对 灰色 模型 的 拟 合 残 差 进行 预报 ， 也 起 到 对 灰色 模型 修 
正 的 作用 ， 实 现 优势 互补 。 
2.2 原始 钟 差 序列 长 度 对 建 模 预 报 结果 的 影响 

不 同 模型 在 预报 时 ， 模 型 参数 主要 依据 最 小 二 乘 原理 解 算 ,因此 建 模 序列 的 长 度 对 模型 
参数 的 计算 会 产生 一 定 的 影响 , 不 同 的 模型 参数 也 会 导致 模型 的 预报 结果 有 所 区 别 。 为 分 析 
建 模 序列 长 度 对 钟 差 预报 精度 的 影响 ， 分 别 以 前 一 天 2h-24h 的 钟 差 序列 构建 二 次 多 项 式 模 
型 、 灰 色 模 型 和 灰色 + 自 回 归 模 型 对 PRN02，PRN06，PRN08 和 PRN29 4 颗 卫 星 接 下 来 
24h 的 钟 差 进行 预报 ， 重 复 预报 30 次 ， 并 以 平均 值 和 均 方 根 误差 作为 精度 评价 指标 。 各 卫 
星 采 用 不 同 长 度 钟 差 序列 预报 的 结果 统计 见 图 6~ 图 9。 
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图 6 PRNO2 不 同 长 度 钟 差 序 列 建 模 预报 精度 统计 


Fig.6 Prediction Accuracy Statistics of Different length Modeling Series of PRNO2 
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图 7 PRN06 不 同 长 度 钟 差 序列 建 模 预报 精度 统计 


Fig.7 Prediction Accuracy Statistics of Different length Modeling Series of PRN06 
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图 8 PRNOS 不 同 长 度 钟 差 序列 建 模 预 报 精度 统计 
Fig.8 Prediction Accuracy Statistics of Different length Modeling Series of PRNO8 
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图 9 PRN29 不 同 长 度 钟 差 序列 建 模 预报 精度 统计 


Fig.9 Prediction Accuracy Statistics of Different length Modeling Series of PRN29 


对 预报 结果 统计 分 析 可 知 : 


(1) 不 同 模型 的 预报 精度 与 原始 建 模 序列 有 一 定 关 系 ， 其 中 二 次 多 项 式 模型 受 其 影响 


最 大 , 当 建 模 序 列 较 短 时 , 预报 误差 较 大 , 随 着 建 模 序 列 增长 , 预报 误差 基本 上 会 逐渐 变 小 ; 


灰色 模型 和 灰色 + 目 回 归 模 型 受 建 模 序 列 长 度 的 影响 相对 较 小 ， 模 型 抗 干扰 能 力 强 ， 预 报 结 


F 
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果 较 为 稳定 ， 这 也 说 明 灰 色 模 型 在 数据 量 较 少 的 情况 下 建 模 效果 比较 好 ， 但 灰色 + 自 回 归 模 
型 的 预报 效果 优 于 灰色 模型 。 以 PRN02 为 例 ， 当 建 模 序列 长 度 为 2h 时 ， 二 次 多 项 式 模 型 和 
灰色 + 自 回归 模型 预报 误差 统计 的 平均 值 和 均 方 根 误差 分 别 是 7.235ns，0.775ns 和 8.743ns、 
0.974ns， 当 建 模 序 列 长 度 为 24h 时 ， 二 次 多 项 式 模型 和 灰色 + 自 回 归 模 型 预报 误差 统计 的 平 
均值 和 均 方 根 误差 分 别 是 1.821ns, 1.093ns 和 1.823ns，1.160ns， 与 2h 的 建 模 序列 预报 结 
相 比 ， 二 次 多 项 式 模型 的 预报 精度 提高 了 79.15%， 灰 色 + 自 回归 模型 的 预报 精度 降低 了 
19.10%， 另 外 ， 对 于 2h 和 24h 的 钟 差 序 列 建 模 预报 结果 ， 灰 色 + 自 回归 模型 的 精度 较 二 次 
多 项 式 模型 分 别 高 了 88.86% 和 36.37%， 说 明 二 次 多 项 式 模 型 适合 原始 数据 量 较 多 时 的 建 模 
预报 ， 而 灰色 + 自 回归 模型 受 原始 建 模 序列 长 度 的 影响 较 小 ， 预 报效 果 较 为 稳定 。 

(2) 不 同 卫 星 预 报 24h 的 最 佳 建 模 序列 长 度 因 使 用 模型 的 不 同 而 不 同 。 对 二 次 多 项 式 
模型 的 预报 结果 分 析 可 知 ,，PRN02 在 序列 长 度 为 4h 时 的 预报 效果 最 佳 ， 而 PRN06，PRN08 
和 PRN29 的 最 佳 序列 长 度 分 别 为 6h，24h 和 22h; 对 于 灰色 模型 ，PRN02，PRN06，PRN08 
和 PRN29 的 最 佳 序列 长 度 分 别 为 4h,2h,1$h 和 24h; 对 于 灰色 + 自 回 归 模 型 ,PRN02,PRN06， 
PRN08 和 PRN29 的 最 佳 序列 长 度 分 别 为 22，2h，15h 和 24h。 

G) 锦 钟 受 建 模 序列 长 度 的 影响 总 体 上 小 于 饮 钟 ， 当 原始 钟 差 序列 长 度 增 加 到 一 定 程 
> E, Huppa té phie E R EET AE o 
le, 2.3 原始 钟 差 序列 波动 对 预报 结果 的 影响 


通过 卫星 钟 差 的 趋势 变化 情况 可 知 , 卫星 钟 差 具 有 绥 慢 下 降 或 者 上 升 的 趋势 , 但 部 分 卫 
© 星 的 钟 差 也 会 在 一 些 时 段 内 波动 较 大 ,对 2016 年 10 月 24 日 的 钟 差 数 据 波动 情况 分 析 可 知 ， 
Ə PRN12 (IIR-M RB) 和 PRN21 OIR RB) 的 钟 差 序列 存在 波动 情况 ， 且 PRNI2 的 波动 幅度 
LO 明显 大 于 PRN21, PAL PÆRER ILE 10， 为 研究 原始 钟 差 序列 波动 对 预报 结果 的 影 
2 响 ， 将 其 与 钟 差 序列 平稳 变化 时 间 段 的 预报 情况 进行 对 比 。 通 过 对 2016 年 10 月 24 日 附近 
N 时 间 段 内 IGS 发 布 的 PRN12 和 PRN21 的 钟 差 产品 比较 , 发 现 2016 年 6 月 27 日 的 钟 差 序列 
I 波动 较为 平稳 ( 见 图 11)， 分 别 采 用 二 次 多 项 式 模型 、 灰 色 模 型 和 灰色 + 上 自 回归 模型 对 10 月 
24 日 和 6 月 27 日 后 一 天 的 PRN12 和 PRN21 的 钟 差 进行 预报 ， 为 便于 对 比 ， 预 报 均 以 24h 
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的 钟 差 序列 建 模 ,图 12 和 图 13 给 出 了 两 种 情况 在 不 同 模型 下 的 预报 误差 ,PRN12; 4I PRN21; 
表示 钟 差 序 列 具有 一 定 波动 ，PRN12, 和 PRN21, 表示 钟 差 序列 变化 趋势 较为 平稳 ， 图 中 从 
上 到 下 分 别 为 二 次 多 项 式 模型 、 灰 色 模 型 和 灰色 + 自 回 归 模 型 的 预报 误差 ， 表 2 给 出 了 两 种 
情况 的 预报 误差 统计 。 
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图 10 PRNI2 和 PRN21 钟 差 序列 (2016 4E 10 H 24 H) 


Fig.10 Satellite Clock Bias of PRN12 and PRN21 (October 24, 2016) 
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图 11 PRNI2 和 PRN21 钟 差 序列 (2016 £ 6 H 27 H) 


Fig.11 Satellite Clock Bias of PRN12 and PRN21 (June 27, 2016) 
对 图 12~ 图 13 和 表 2 分 析 可 知 : 
(1) 原始 钟 差 序列 波动 对 预报 结果 有 一 定 影响 ， 对 比 不 同 模型 的 预报 结果 可 以 看 出 ， 
当 钟 差 序 列 变化 趋势 平稳 时 ， 预 报 误 差 总 体 上 变化 较为 平缓 ， 且 变化 幅度 小 ， 预 报 精度 高 。 
PRN12; 和 PRN12; 采用 二 次 多 项 式 模 型 、 灰 色 模 型 和 灰色 + 自 回归 模型 预报 误差 的 RMS 均 
方 根 误差 PRN12 分 别提 高 了 18.6%，32.0% 和 9.6%; PRN21; fI PRN21; 采用 二 次 多 项 式 模 
型 、 灰 色 模 型 和 灰色 + 自 回 归 模 型 预报 误差 的 均 方 根 误差 分 别 是 2.191ns，2.572ns，2.572ns 
和 1.675ns, 1.261ns, 1.240ns, PRN21; 的 预报 精度 较 PRN21; 分 别提 高 了 23.696, 51.0968 
51.896. 
(2) 不 同 模型 的 预报 结果 与 钟 差 序列 的 波动 幅度 大 小 有 关 。 当 钟 差 序列 波动 较 大 时 ， 
灰色 + 自 回 归 模 型 的 预报 效果 较 好 ，PRN121 采用 灰色 + 自 回 归 模 型 的 预报 精度 较 二 次 多 项 式 
模型 和 灰色 模型 分 别 高 了 36.290651 34.9%， 当 钟 差 序列 波动 幅度 不 大 时 ， 二 次 多 项 式 模型 的 
预报 效果 较 好 , PRN21, 采用 二 次 多 项 式 模 型 的 预报 精度 较 灰 色 模 型 和 灰色 + 自 回 归 模 型 分 别 
高 了 14.8% 和 14.896. 
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de 2 钟 差 波 动 时 预报 结果 统计 《〈 单 位: ns) 


Table 2 Statistics of Prediction Results when Clock Bias Fluctuates (unit:ns) 


MAX MIN MEAN RMS 

QP 0.588 -1.862 -0.408 0.716 

PRNI2, GM(1,1) 1.140 -1.561 -0.113 0.702 
GM(1,1)-AR 1.140 -1.560 -0.085 0.457 

QP 0.524 -0.913 -0.281 0.583 

PRN122> GM(1,1) 0.173 -1.281 -0.633 0.477 
GM(1,D)-AR 0.170 -1.282 -0.652 0.413 


PRN21, QP 4.083 -0.535 1.749 2.191 


GM(1,1) 4.754 -0.418 2.098 2.572 


GM(1,I)-AR 4.754 -0.420 2.097 2.572 

QP 1.128 -3.913 -1.114 1.675 

PRN215 GM(1,1) 2.161 -0.047 1.182 1.261 
GM(1,1)-AR 2.161 -0.556 1.146 1.240 


3 结论 


本 文选 取 3 种 预报 模型 从 不 同 角度 对 全 球 定位 系统 星 载 原 子 钟 的 钟 差 预报 进行 了 研究 ， 
结果 表明 : 

OD 从 不 同类 型 原子 钟 的 预报 结果 可 知 ， 铭 钟 短期 预报 的 稳定 性 总 体 上 优 于 钨 钟 ， 但 
根据 预报 模型 不 同 ， 钨 钟 的 预报 精度 也 可 与 鲍 钟 相当 。 灰 色 模 型 和 灰色 + 自 回 归 模 型 对 于 锦 
钟 和 钨 钟 的 预报 结果 均 比较 稳定 ， 误 差 变 化 幅度 不 大 ， 且 预报 精度 相差 不 多 , 但 二 次 多 项 式 
模型 对 钨 钟 的 预报 误差 波动 幅度 较 大 , 而 锦 钟 预报 误差 的 波动 幅度 较 小 , 可 与 灰色 模型 和 灰 
色 + 自 回归 模型 相当 。 总 体 来 说 ， 当 建 模 序 列 为 24h 时 ， 采 用 灰色 + 自 回 归 模 型 的 预报 效果 
最 好 。 


(2) 不 同 模型 、 不 同 长 度 的 钟 差 序 列 建 模 预 报效 果 不 同 ， 二 次 多 项 式 模型 受 建 模 序 列 
长 度 的 影响 最 大 ， 灰色 + 自 回归 模型 受到 了 影响 最 小 ， 灰 色 模 型 的 预报 效果 随 建 模 序列 长 度 
变化 的 趋势 与 灰色 + 自问 归 模 型 相似 ;采取 不 同 模型 预报 时 ， 不 同类 型 原子 钟 的 最 佳 序列 长 
度 也 有 所 区 别 ; 另外 , 在 一 定 程 度 下 ， 钢 钟 受 建 模 序 列 长 度 的 影响 小 于 钨 钟 ， 但 当 建 模 序 列 
增加 到 一 定 程度 时 ， 锦 钟 和 钨 钟 的 预报 精度 会 趋 于 稳定 。 

(3) 当 建 模 序 列 长 度 为 24h 时 ， 钟 差 序 列 变化 趋势 稳定 时 的 预报 结果 优 于 钟 差 序列 波 
动 时 的 结果 ， 当 钟 差 波动 幅度 不 同时 ， 不 同 预 报 模型 受到 的 影响 也 有 所 区 别 。 
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Research on the influence of prediction model and 
modeling sequence length on short-term prediction 


accuracy of clock bias 


GUO Zhongchen!, SUN Peng, LI Zhichun', BAI Hongwei' 
(1. School of Environment and Surveying Engineering, Suzhou University, Suzhou, Anhui 234000, China ) 

Abstract: Aiming at the problem that different types of GPS space-borne atomic clockshave 
different adaptability to different prediction models in clock difference prediction, QP model, 
GM(1,1) model and GM(1,1)--AR model are used to predict the clock bias of different types of 
atomic clocks. The analysis focuses on the prediction accuracy of different types of atomic clocks, 
the effect of clock bias series of different lengths on model modeling and forecast, and the effect 
of clock sequence fluctuations on model forecasting. The results show that the accuracy of the 
clock bias prediction has a certain relationship with the length of the modeling series. The QP 
model is most affected, and the GM(1,1)+AR model is least affected. The optimal modeling 
sequence length for different types of atomic clocks is different under different prediction models. 
The rubidium clock is less affected by the length of the modeling series than the cesium clock. QP 
model has poor prediction effect on cesium clock, and the prediction effect on rubidium clock is 
comparable to GM (1,1) model and GM (1,1) + AR model. When the clock bias series fluctuates, 
the accuracy of modeling prediction decreases, and the prediction results of different models are 
affected differently by the amplitude of the clock difference fluctuation. 

Keywords: clock bias; short-term prediction; modeling series length; clock bias fluctuation; 


Prediction model 
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